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Compression index (Cc) of normally consolidated (NC) clays 
determined by the oedometer experiments is utilized for calculating 
the consolidation settlement of shallow foundations. The 
determination of the Cc from the tests takes much more time and 
so empirical correlations based on clay properties can be a suitable 
alternative for the prediction of settlement. However, uncertainty in 
the measurements of input parameters has always been a major 
concern. Robust optimization is very popular due to its 
computational tractability for many classes of uncertainty sets and 
problem types. Therefore, in this research, an innovative method 
based on robust optimization has been used to investigate the effect 
of such uncertainties. To achieve these, the results of 433 
oedometer tests taken from geotechnical investigation locations in 
Mazandaran province of Iran have been used. Based on Frobenius 
norm of the data points, uncertainty definition is presented and 
examined against the correlation coefficients for several empirical 
models for predicting Cc value and thus optimum values are 
determined. The results in compare with previous models indicate 
the robust method is a better pattern recognition tool for datasets 
with degrees of uncertainty. The variation of the Cc values with 
soil properties, namely, water content (w,,), initial void ratio (e,), 
and liquid limit (ZL), by considering different value of 
uncertainties (0, 5 and 10%) was considered and indicated that the 
effect of e, is more than other two physical parameters (w, and 
LL). The best model performance during in deterministic valuation 
and considering uncertainty is further shown. 
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1. Introduction 


Stress increases in soil layers owing to structures are accompanied by settlements. Settlement- 
based design is an essential and important criterion for shallow foundations [1,2] The volume 
change characteristics of clays can be defined by oedometer tests using Terzaghi's one- 
dimensional consolidation theory [3]. Compression index (C,), a compressibility parameter, 
presents the slope of the curve of the void ratio (e) versus the logarithm of the effective pressure 
and is obtained by these tests. Additionally, the C. value, the main indicator of the 
compressibility of soil, is generally utilized for calculating the settlement of soil layers [1,4]. C. 
is an important parameter used for the estimation of the primary consolidation settlement of clays 
[5]. High value of C;, indicates large settlement. The calculation of C, value from oedemeter tests 
is relatively time consuming and expensive [6-9]. Performing oedemeter tests also requires 
precision, precautions and expertise and so is cumbersome [5]. Even a small disturbance can 
cause overestimation or underestimation of the C. value and so the settlement [5]. Thus, if 
correlations between the C. value and easily determined soil properties are developed, the 
determination of the C. value may save a lot of time and laboratory expenditure [5]. Therefore, 
many investigators have performed various investigations to predict the C, value from the soil 
properties determined more easily. They suggested single or multi-variable equations for the 
prediction of the C, value, which is presented in Table 1. In this table, LL is the liquid limit; w,, is 
natural water content; e, is the in-situ void ratio; G, is the specific gravity, and yq is the dry unit 
weight. All estimations based on the empirical methods need assessment of input parameters 
which incline to uncertainties and inaccuracies. Thus, a formula that is capable of resolving the 
uncertainties and inaccuracies in input parameters requires to be suggested. 


Data uncertainty is observed generally in many optimization problems [10]. Measurement errors 
can lead to uncertainty in the parameters of an optimization model. Recently, robust optimization 
(RO) approach, a methodology to solve mathematical optimization problems with uncertain data, 
is developed to address such uncertainties [11] The objective of the RO is to find solutions that 
are immune to all perturbations of the data in a so-called uncertainty set. RO is popular because 
it is a computationally tractable methodology and has a wide range of applications in practice 
[12]. 


Keeping this in view, a RO approach was proposed in this paper in order to quantify the effect of 
uncertainties on evaluation of the input parameters. The RO model is the robust counterpart of 
the least square model and is a second order cone program (SOCP) in which, possible 
uncertainties can be reasonably adjusted. SOCP has widely been used in operation research [13]. 
In geotechnical engineering field, Kalantary et al. [14] and MolaAbasi et al. [15] developed 
robust optimization model to predict the liquefaction induced lateral displacement and shear 
wave velocity, respectively. It is shown that if uncertainty (%) is set to zero, the method reduces 
to ordinary regression method. Also, Kalantary et al. [14] have shown that logarithmic lateral 
spreading correlations perform better in deterministic valuation, whereas by considering 
uncertainty, they give similar degrees of accuracy to new linear model. Thus, this method can be 
used in empirical correlation of C, as a step forward in comparison with statistical approaches 
that considers the variation of uncertainty. 
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Table 1 
Some widely used compression index equations. 
Indepe gent Equation References Fauaon 
variable no. 
C, = 0.010, — 0.05 Azzouz et al. [16] (1) 
es C, = 0.010, Koppula [17] (2) 
me C. = 0.01w,, — 0.075 Herrero [18] (3) 
be C. = 0.0130, — 0.115 Park and Lee [19] (4) 
& C, = 0.49) — 0.11 Park and Lee [19] (5) 
g C, = 0.4(e9 — 0.25) Azzouz et al. [16] (6) 
3 ‘i C, = 0.287e,) — 0.015 Ahadiyan et al. [20] (7) 
oO 0 
ra C, = 1.02 — 0.95e, meme Arman (8) 
& 
S$ C, = 0.287e, — 0.015 Farzi [22] (9) 
= C, = 0.006(LL — 9) Azzouz et al. [16] (10) 
5 C. = (LL — 13)/109 Mayne [23] (11) 
LL C, = 0.009(LL — 10) eae ts Peck (12) 
C, = 0.014LL — 0.168 Park and Lee [19] (13) 
PI C, = 0.0082PI + 0.0475 Jain et al. [25] (14) 
LL 
LL,Gs ve= 0.2926(—~) 6, Park and Lee [19] (15) 
wo... LL C. = 0.009w, + 0.005LL Koppula [17] (16) 
2 a C, = 0.009w, + 0.002LL — 0.1 Azzouz et al. [16] (17) 
5 C0, Wn C. = 0.4(€9 + 0.001, — 0.25) Azzouz et al. [16] (18) 
5 Al-Khafaji and 
s Pe Gs C. = —0.156 + 0.411e9 + 0.00058LL Andersland [26] (19) 
5 C, = —0.023 + 0.271e, + 0.001L Ahadiyan et al. [20] (20) 
. €05 Wns C. = 0.37(€9 + 0.003LL + 0.0004w,, — 0.34) Azzouz et al. [16] (21) 
+ LL C. = —0.404 + 0.341e, + 0.006w, + 0.004LL Yoon and Kim [27] (22) 
E Goo  C, = 0.141G24[(1 + e9)/G.]?38 Herrero [18] (23) 
C,. = 0:1597 (a, 9" 4771 
= = Ozer et al. [28 24 
Wns + €y)1592 (LL 0.0638) (770-8276) [28] (24) 
LL,€9,Ya  C, = 0.151 + 0.001225w,, + 0.193e9 Ozer et al. [28] (25) 


— 0.000258LL — 0.0699y, 


The study presented here aims both to predict the Cc values from the easily determined soil 
parameters by using RO model without performing consolidation tests and to evaluate the effects 
of uncertainty of each parameter (wn, eo, Gs and LL) independently on the analysis outcome Cc. 
To achieve these, the results of 433 laboratory tests (oedometer) on clay at 115 projects in north 
of Iran (Mazandaran province) were used while developing RO model. The measured C, values 
from the experiments were compared with the predicted C. values with different level of 
uncertainties. 


The rest of the paper is arranged as follows. A review of RO model is given in Section 2. While 
the database compilation is presented in Section 3, the development of the ten different RO 
models is given in Section 4. The results of ten RO models and the discussion the RO results 
with the previous models in the literature given in Section 5. 
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2. Review of robust optimization (RO) model 


In mathematical optimization models, it is generally assumed that the input data is precise and 
the effect of uncertainties on the optimality and feasibility of the models are ignored. It is thus 
possible that when the data is different from the assumed nominal values, the generated optimal 
solution may violate critical constraints and make it poor from an objective function point of 
view. These explanations provoke the need for methodologies in mathematical optimization 
models that explain solutions immune to data uncertainty [29,30]. For instance, inaccuracies 
enter in field measurements of C. in the case histories alike all other natural phenomenon 
measurements. Such inaccuracies exist in other influencing parameters and can cause deviation. 
If such deviations are given as boundaries of the central point of the data (Figure 1) and the true 
data point could exist at any point within this boundary, the robust optimization (RO) 
investigates to minimize the maximum error with any particular level of uncertainty. In the 
particular case of zero uncertainty, RO approach decreases to normal regression analysis. 


. 
Robust 
Method 
© t 
=} 
oO 
=> 
— . 
= 
2 
= 
oO 
Regression Boundary 
Method 
e Actual Datapoint 


Input Value 
Fig. 1. Sample view of uncertainties. 


In the past decade, there were noticeable developments in the theory of robust convex 
optimization. In this paper, a RO model for the least squares method was considered. It was 
firstly assumed that the level of uncertainty of a database is unknown and equal to p. Then, the 
RO model takes into account this level of uncertainty in the database and minimizes the worst 
case residual 1.e. 


MIN, MAX ¢ -1) <p || (A+ E)x-(b+r)|| (26) 
Where E and r are uncertainties in A and b, respectively, and the matrix norm, || || is the 
Frobenius norm. Definitely, Eq. (26) cannot be solved by utilizing classical optimization 


algorithms. However, it can be rewritten in the following second order cone program (SOCP) 
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form [31] (Eq. 26) and solved using an effective software such as SeDuMi [26], which is an 
interior point based software for solving SOCP and semi definite optimization. 


min(t + ps) (27) 


with the constrains, ||Ax — b|| < t, and /1 + ||x||? < s 


First we rewrite problem (26) in the dual form of SeDuMi’s input format namely 


max b’y 
c-AT’y EK (28) 
where 
0 -1 0 Oixn 
—b Omx1 = Omx1 -A -1 t 
c-{ { |,at=[ 0 -1 Onn |, B=(?), ¥=(9) (29) 
Onx1 0 0 Oishi nxt 
Onx1 Onx1 —Inxn 
K=O XOua (30) 


and Q;, denotes the second order cone in R*. 


This kind of SOCP problem can be easily solved by means of SeDuMi which is one the most 
efficient interior-point methods based software packages for solving SOCP for four different 
values of uncertainty parameter. Thereafter, SeDuMi called by the following MATLAB 
command was used to solve it four different values of p in order to evaluate uncertainty: 


[x,y|=sedumi(At,b,c,K) 


Where At denotes transpose of the matrix A® in Eq. (29); b and c are taken from Eq. (29); and K 
also is obtained from Eq. (30). The x and y outputs denote the solutions of Eq. (30) and its dual 
problem. Then, uncertainty is calculated from the following equation: 


Uncertanity (%) = TaN PS x 100 (31) 


where ||DATA||;-. is the Frobenius norm of maximum error of data matrix. 


3. Database compilation 


Databases have compiled the data from 433 oedometer tests for clays at 115 building sites in 
north of Iran (Mazandaran province), which is taken from the Geotechnical Report of Site 
investigations performed at Technical and Soil Laboratory of Mazandaran Province between 
2000-2010 years. Considering the past investigations in the literature, in this paper, the C, value 
of the clays was expected to be influences by the initial void ratio (e,), liquid limit (LZ), natural 
water content (w,,), plastic index (PJ), and specific gravity (G;). The locations of collected 


6 Y. Erzin et al./ Journal of Soft Computing in Civil Engineering 4-3 (2020) 01-16 


database contain 433 data summarized in Figure 2. The samples were all collected and tested 
according to standard ASTM D2435-96. The database covers a wide spectrum of fine grained 
soils and soil parameters, including e,, w,, LL, G,, and C,. The details of soil parameters used 
while developing the RO models are given in Figure 3 and in Table 2. 


Caspian sea 


Golestan 


Fig. 2. Descriptive data collection location. 
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0 : 0 0.0 2.40 0.0 : hae 
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Fig. 3. Descriptive statistics of variables. 
Table 2 
The details of the parameters used in the RO models developed 
Parameters used Minimum Maximum Mean Std. Deviation 
Input parameters 
w (%) 5.00 70.00 28.6395 7.8835 
Ya (kKN/m?) 9.30 22.4 15.048 1.3942 
LL (%) 24 81 39.8819 9.7299 
G, 2.43 2.80 2.6385 0.0603 
e 0.18 1.88 0.7693 0.1796 


Output parameter 
C. 0.05 0.63 0.2084 0.0823 
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4. Development of robust optimization (RO) models 


In order to investigate the effect of uncertainties in predicting the Cc, value ten different Robust 
Optimization (RO) models were developed by considering the previous studies given in Table 1. 
Description of the RO models used in this study is given in Table 3. Evaluating a; coefficients in 
the ten different RO models given in Table 4 by means of robust optimization method constitutes 
the main goal of this paper. 


Table 3 
Description of different RO models used in this study. 
Model Model description 
ROI Co = AyWn + Az 
RO2 Co = ay,€y + Az 
RO3 C, = a,LL +a, 
RO4 C, = a,LL™’Gs® 
ROS Co = Ay Wn + AQZLL + a3 
RO6 C. = a,€y + aAZLL + a3 
RO7 Co = Ayn + AQLL + A3Wy + A4 
RO8 Co = 4G" [(1 + €9)/Gs] 8 
ROS Co = ay (wn?) (1 + €0) LL"*)(¥q°) 
RO10 C. = A, + A2W, + AZLL + Ayeg + AsVq 
Table 4 
The determined coefficients for ten different RO models. 
Model ay a a a4 a5 
ROI 0.0931 0.1700 - - - 
RO2 0.146 0.0904 - - - 
RO3 0.0879 0.161 - - - 
RO4 0.65726368 0.40297 -0.40743 - - 
ROS 0.0848 0.0805 0.143 - - 
RO6 0.133 0.0506 0.0794 - - 
RO7 0.125 0.0479 0.0479 0.0744 - 
RO8 0.621015 -0.46271 0.21201 - - 
RO9 0.714573 0.46383 -0.16917 0.32257 -0.15777 


RO10 0.051894 0.036503 0.035216 0.094971 0.034755 
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Therefore, these ten RO models can be formulated of the following form 
Ax =b,A€R™", be R™1 (32) 


which shows over-determined set of equations in the case of (m>n). Initially, all the RO models 
are examined using the compiled database and simple regression analyses, containing least 
square approaches. The classical method for solving least squares problems, e.g.: 


Minyepnr||Ax — b|| (33) 


In order to evaluate the prediction performance of the RO models developed in this study, 
different performance indices were used [32—34]. Four performance indices, namely, absolute 
fraction of variance (R?), given by Eq. (34), root mean squared error (RMSE), given by Eq. (35), 
mean absolute percent error (WAPE), given by Eq. (36), and mean absolute deviation (MAD), 
given by Eq.(37), were calculated for examining of the prediction capacity of each RO model. 


on Be “ 
RMSE = EM Cn — Gi)? (35) 
MAPE = al x 100 (36) 
MAD = zien (37) 


where / is the number of data; and C,,; and C,; are observed and calculated values, respectively. 
The better model performance will be obtained by lower RMSE, MAPE and MAD values. Under 
ideal conditions an accurate and precise method gives R* of 1.0, RMSE, MAPE, and MAD of 
zero. 


5. Results and discussions 


As mentioned earlier, to evaluate the uncertainty, variations of the coefficients (a), a2, a3, a4, and 
as) presented in Table 3 have been evaluated against various levels of uncertainty for ten RO 
models. The results are summarized in Figures 4 and 5.These figures show that by increasing the 
uncertainty, the variability approaches a nearly stable value as observed by Kalantary et al. [14]. 
It must be pointed out at this stage that if uncertainty (%) is set to zero, the method reduces to an 
ordinary multiple linear regression (MLR) technique. The coefficients determined from Figs. 4 
and 5 are given in Table 4 for ten RO models. As known, the physical parameters of soils have 
an important influence on soil’s compressibility parameters. The variation of the C, values with 
soil properties (w,, e,, and LL) by considering different value of uncertainties 0, 5 and 10% is 
also given in Figure 6. When compared a, coefficients given in Table 4 for the ROI, RO2 and 
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RO3 models developed using single variable, the RO3 model yielded the highest a; coefficient. 
This result indicates that the effect of e) is more than other two physical parameters (@, and LL). 
The effect of e9 was also observed in the models RO6, RO7 and R10 developed using two, three 
and four variables, respectively. 


In order to determine the precision of each RO model, the performance indices (R’, RMSE, 
MAPE, and MAD) calculated by using Eqs. (34) to (37), respectively, were calculated for each 
RO model with different value of uncertainties, which is shown in Figure 7. As mentioned 
earlier, if uncertainty (%) is set to zero, the method reduces to an ordinary multiple linear 
regressions (MLR) technique. The calculated performance indices (R’, RMSE, MAPE, and MAD) 
for the uncertainty of 100% are also given in Table 5 for ten RO models. The models RO2, RO6, 
and ROT yielded lower RMSE, MAPE, and MAD values and higher R’ values than the other RO 
models (Table 5). 


0.8 7 0.4 
0.7 =e al a2 0.35 1 al a2 
0.6 0.3 
0.25 
z 0.5 2 ai 
504 Ss 
i] = 0.15 
03 5 
i) S 0.1 
{S) 02 1S) 
0.05 
ai 4 | 0 
0 -0.05 
Model RO1 Model RO2 
0.1 -0.1 
0 20 40 60 80 100 0 20 40 60 80 100 
Uncertainty (%) Uncertainty (%) 
0.35 1 
- ental | a2 08 | —e—al a2 tr 23 
0.25 04 | : 
= = 02 
g % Pia 
y i) 0 i } } 
Bo.15 5 ‘i | 
'S) oO 
a | - es es er en | 
-0.6 += os T 
0.05 ’ f 
Model RO3 28 7 | Model RO4 
0 -] 
0 20 40 60 80 100 0 20 40 60 80 100 
Uncertainty (%) Uncertainty (%) 


Fig. 4. The variation of constant coefficients (a;) versus uncertainties for ROI, RO2, RO3, and RO4 
models. 
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al a2 


a3 


Model ROS 


0 20 40 60 
Uncertainty (%) 


80 100 


Coefficients 


Coefficients 


Coefficients 


Uncertainty (%) 


0.35 
0.3 
0.25 
0.2 
0.15 
0.1 


0.4 ———al a2 —t— 93 
0 
0 20 40 60 80 100 
Uncertainty (%) 


Coefficients 


Uncertainty (%) 


Fig. 5. The variation of constant coefficients (a;) versus uncertainties for ROS, RO6, RO7, RO8, ROS, 
and RO10 models. 
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Fig. 6. The variation of C, outputs with soil properties (w,, e,, and LL) by considering different value of 


uncertainties of 0, 5 and 20%. 
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Fig. 7. The performance indices (R’, RMSE, MAPE and MAD) versus uncertainties in predicting C, value 


Table 5 


for different ten RO models. 


The performance indices calculated for all RO models. 

Model R? RMSE MAPE MAD 
RO1 0.8788 0.0780 26.9115 0.0561 
RO2 0.9191 0.0637 22.6891 0.0473 
RO3 0.8717 0.0803 27.7060 0.0577 
RO4 0.7019 0.1223 51.9701 0.1083 
RO5 0.8865 0.0755 26.3032 0.0548 
RO6 0.9194 0.0636 22.5631 0.0470 
RO7 0.9226 0.0623 22.2321 0.0463 
RO8 0.4031 0.1731 77.4982 0.1615 
RO9 0.8900 0.0743 29.8022 0.0621 
RO10 0.9019 0.0702 24.7954 0.0517 


In order to compare the performance of RO models with the correlations suggested by previous 
researchers and given in Table 1, the performance indices (R’, RMSE, MAPE, and MAD) were 
calculated for these correlations by using the test data in this study which is depicted in Table 5. 
When compared the performance indices of the correlations developed by previous researchers 
using single variable of @, (i.e. Eqs. 1 to 4) with those of the ROI model developed by using the 
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same variable, the correlations developed by Azzouz et al. [16] and Herrero et al. [18] 
represented by Eqs. (2) and (4), respectively, yielded higher R? values and lower RMSE, MAPE, 
and MAD values. When compared the performance indices of the correlations developed by 
previous researchers using single variable of eg (i.e. Eqs. 5 to 8) with those of the RO2 model 
developed by using the same variable, higher R? values and lower RMSE, MAPE, and MAD 
values were obtained for the correlations developed by Azzouz et al. [16] and Ahadiyan et al. 
[20], represented by Eqs. (6) and (7), respectively, and the RO2 model. When the performance 
indices of the correlations developed by previous researchers using single variable of LL (i.e. 
Eqs. 10 to 13) were compared with those of the RO3 model developed by using the same 
variable, the correlation developed by Azzouz et al. [16], represented by Eq. (10), and the RO3 
model produced higher R’ values and lower RMSE, MAPE, and MAD values. It can be noted 
from Tables 5 and 6 that the highest performance was obtained for the correlation suggested by 
Azzouz et al. [16], (Eq. 6), by using single variable of eo, which indicates the higher effect of e9 
on the compression index (C, ) value as mentioned earlier. 


Table 6 
The performance indices calculated for the best used compression index equations. 
Independent variable References Equationno R? RMSE MAPE MAD 
& Azzouz et al. [16] (1) 0.9105 0.0670 26.0071 0.0542 
a= ii 
= Herrero [18] (3) 0.9259 0.0610 22.6792 0.0473 
oO 
oO 
S Azzouz et al. [16] (6) 0.9490 0.0506 18.7712 0.0391 
= = 
= Ahadiyan et al. [20] (7) 0.9461 0.0136 19.0337 0.0397 
“oh 
5 LL Azzouz et al. [16] (10) 0.8577 0.0845 30.8510 0.0643 
Wy» LL Azzouz et al. [16] (17) 0.9141 0.0657 25.6741 0.0535 

mn £0, Wn Azzouz et al. [16] (18) 0.9460 0.0521 19.5910 0.0408 
= 
z Al-Khafaji and Andersland [26] (19) 0.9060 0.0687 26.1167 0.0544 
= Co, LL 
5 Ahadiyan et al. [20] (20) 0.9344 0.0574 19.9235 0.0415 
ia) 
= Co, Wn, LL Azzouz et al. [16] (21) 0.9520 0.0491 18.0270 0.0376 
ma 
= G,, eo Herrero [18] (23) 0.9274 0.0604 21.0635 0.0439 
= 

Ozer et al. [28] (24) 0.9482 0.0510 19.3152 0.0402 

@,, LL, €0, Ya 
Ozer et al. [28] (25) 0.9441 0.0530 20.2145 0.0421 


When compared the performance indices of the correlation developed by Park and Lee [19](Eq. 
4) using two variables of LL and G, with those of the RO4 model developed by using the same 
variables, both the correlation and the RO4 model yielded lower R’ values and higher RMSE, 
MAPE, and MAD values, which was not acceptable. When the performance indices of the 
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correlations developed by past researchers using two variables of w, and LL, Eqs. (16) and (17), 
were contrasted with those of the ROS model developed by using the same variables, the 
correlation developed by Azzouz et al. [16], represented by Eq. (17), yielded higher R’ values 
and lower RMSE, MAPE, and MAD values. When compared the performance indices of the 
correlations developed by past researchers using two variables of eg and LL, Eqs. (19) and (20), 
with those of the RO6 model developed by using the same variables, the correlation developed 
by Ahadiyan et al. [20], represented by Eq. (20), yielded higher R? values and lower RMSE, 
MAPE, and MAD values. When compared the performance indices of the correlation developed 
by Herrero et al. [18] using two variables of e9 and G;, Eq. (23), with those of the RO8 model 
developed by using the same variables, both the correlation yielded higher R? values and lower 
RMSE, MAPE, and MAD values. It can be noted from Tables 5 and 6 that when compared the 
prediction performance of the correlations and RO models developed by using two variables, the 
highest performance obtained for the correlation suggested by Azzouz et al. [16], (Eq. 18), by 
using two variables of eg and w,, which indicates the higher effect of e9 on the compression index 
(C.) value as mentioned earlier. 


When compared the performance indices of the correlations developed by past researchers using 
four variables of eo, @n, ya, and LL, Eqs. (24) and (25), with those of the models RO9 and RO10 
developed by using the same variables, the model RO10 and both correlations developed by 
Ozer et al. [28] yielded higher R’ values and lower RMSE, MAPE, and MAD values. 


When compared the performance of ten different Robust Optimization models, developed for the 
examination of the influences of uncertainty of each soil parameter independently on the analysis 
outcome Cc, the robust method is a better pattern recognition tool for datasets with degrees of 
uncertainty. The linear model contains LL (RO3) are found to be more effective than the 
available empirical formulas. The variation of the C. values with soil properties (w,, e@,, and LL) 
by considering different value of uncertainties (0, 5 and 10%) was evaluated and indicated that 
the effect of e is found to be more than other two physical parameters (@, and LL). 


6. Conclusions 


The determination of the Compression index (Cc) of normally consolidated clays utilized to 
assess the consolidation settlement of shallow foundations takes much more time. Therefore, 
empirical correlations between Cc value and soil properties can be a suitable alternative for the 
prediction of Cc value. However, uncertainty in the measurements of input parameters has 
always been a major concern. Robust optimization is widely used due to its computational 
tractability for many classes of uncertainty sets and problem types. Therefore, in this research, an 
innovative method based on robust optimization has been used to investigate the effect of such 
uncertainties in the prediction of the Cc value. To achieve this, 433 data points have been taken 
from geotechnical investigation locations in Mazandaran province of Iran. The identification 
procedure used is based on the robust model of the least square which is a second order cone 
problem and is capably solved by the interior point technique. Based on Frobenius norm of the 
data points, uncertainty definition is presented and examined against the correlation coefficients 
for several empirical models and thus optimum values are determined. The comparison of results 
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of the RO models with those of previous models indicate the robust method is a better pattern 
recognition tool for datasets with degrees of uncertainty. The variation of the Cc values with soil 
properties, namely, water content (w,,), initial void ratio (e,), and liquid limit (ZL), by 
considering different value of uncertainties (0, 5 and 10%) was considered and indicated that the 
effect of e, is more than other two physical parameters (w, and LL). Also, the linear robust 
model developed by using LL values is more efficient than the available empirical formulas. The 
best model performance during in deterministic valuation and considering uncertainty is further 
shown. 


It should be noted that the robust models to predict Cc were obtained using limited amount of 
database. Again, more studies are required to check the validity of the equations that have been 
derived for other database. 
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